Energy concentration and Sommerfeld condition for Helmholtz 
equation with variable index at infinity 
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Abstract 

We consider the Helmholtz equation with a variable index of refraction n(x), which is not 
necessarily constant at infinity but can have an angular dependency like n(x) n^x/lxl) as 
| a; | — > oo. Under some appropriate assumptions on this convergence and on we prove that 
the Sommerfeld condition at infinity still holds true under the explicit form 

dx 



Vu 



inH 2 u 



x 



< +oo. 



Li- 



lt is a very striking and unexpected feature that the index appears in this formula and not 
the gradient of the phase as established by Saito in j2S| and broadly used numerically. This 
apparent contradiction is clarified by the existence of some extra estimates on the energy decay. 
In particular we prove that 
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In fact our main contribution is to show that this can be interpreted as a concentration of the 
energy along the critical lines of n^. In other words, the Sommerfeld condition hides the main 
physical effect arising for a variable n at infinity; energy concentration on lines rather than 
dispersion in all directions. 



Introduction 



We consider the Helmholtz equation with a variable index of refraction n(x), with a slow, and 
only radial decay to a constant n^x/lxl) at infinity 

ieu e + Au £ + n(x)u £ = — f(x), e > 0. (1.1) 

Our main interest is the so called limiting absorption principle (i.e. to study the limit when 
e > approaches to in (jl.lj) ) and the validity of the Sommerfeld radiation condition at 
infinity. One of the main results in this paper is to prove that 
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A direct consequence of this condition is the more classical setting 



liminf 



|x|=r 



Vu(x) — inll 2 (- — :)u(x) -j— r 



"oc 



da(x) — > 0, as r^oo, (1.3) 



where da denotes the Lebesgue measure on the sphere. It is a very striking and unexpected 
feature that the term 

appears in this formula instead of V(p(x) with (p(x) the solution to the corresponding eikonal 
equation 

|Vyf = n, 

as established by Saito in [2H| . This phenomenon, as well as the proof of (jl.3j) . can be explained 
by the existence of some new energy estimate that we state later on and which is, in some sense, 
the main result of this paper. It explains that the Sommerfeld radiation condition hides the 
main physical effect of a variable n at infinity; energy is not dispersed in all directions but 
concentrated on those given by the critical points of Vn(-pj). It would be interesting to prove 

\X\ 

that only local maxima of n can carry energy. 

These estimates use in a strong way the inequality obtained in |2I] for the tangential part 
of the gradient of the solution of In order to recall that result we need some notation. 

Firstly we define, for j G Z, the annulus C(j) by 

C(j) = {x G R d s.t. 2 j < \x\ < 2 J+1 }. 

Then we set 

INIIflo := su p o / H 2 dx i ( 1A ) 

R>R K JB{R) 

■■= J2i 2J+1 1 \f\ 2 dx ^' 2 + / m 2 d A l, \ (1.5) 

with J defined by 2 J < R < 2 J+1 , and we drop the index R if Rq = 0. We also denote the 
radial and tangential derivatives by 

d x d x d 

—u(x) := — • Vu(x), V T u(x) = —u(x) := Vu(x) — — —u, (1.6) 
or \x\ or \x\ or 

and for a function n(u) G C 1 (5 o!_1 ), we shall consider 

v u ,n{w) = —n{u):=\x\—n{ T1 ), uj = — . 

ooo or \x\ \x\ 

Let us consider the following assumptions: 

n = n 1 + n 2 with n 2 G L°°, n > 0, (1.7) 
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|| {ni) l l 2 u\\2 < (1 — co)||Vu|| 2 for all smooth functions u and some Co > 0, (1.8) 

2> sup^ \" :=(3< 1. (1.9) 

S c W "GO 

Above (a)_ denotes the negative part of a 6 R. 
In [21] we prove the following result. 

Theorem 1.1 We assume one of the following two conditions: 

(i) d>3, fnP- fOp and i? = 0; 

^ d = 2, fli.7| )-| fT~^) . n > n > and i? = n^ 1 ^ ; 

Then the solution to the Helmholtz equation hl.l)) satisfies, 

M := IHVulll^, + llnaut + / 1 1 dx 

J\x\>Ro Fl 

<C( £ + ||n 2 || oo )iV ii0 (^) 2 . (1.10) 



The homogeneity of the above estimate makes it compatible with the high frequencies (replace 
n by fi 2 n). Our main interest in [21] was to obtain estimates with the right scaling. In particular 
we were able to recover the well known inequality for u of Agmon and Hormander in [3j in the 
constant coefficient case. Similar results but not scaling invariant were obtained in ^3] and 
. The scaling plays a fundamental role in the applications to nonlinear Schrodinger equations 
(|16j) and in the high frequency limit for Helmholtz equations (jSj, |Zj)- m this paper we get 
an explicit Sommerfeld radiation condition for solutions obtained from the limiting absorption 
principle. As we have already said the estimate of the tangential component of the gradient 



— dx < oo, 1.11, 



\x\>R 



\X\ 



given in (jl.lUJl turns out to be fundamental. In order to get it we need the smallness assumption 
given in ()1.9)) . We do not know if it is necessary or not. However the condition in (J1.9)) is 
necessary and can not be relaxed to a Coulomb type of decay, even if smallness is added. This 
is proved in the appendix using as counterexamples a family of wave guides for which the 
estimate of the tangential derivative is false. 

In order to prove the new energy estimate, we need to impose some extra assumptions on n. 
They are the following ones: 

there exists ri 0O (^) G C 3 (^ d_1 ), noo(A) > n > 0, (1.12) 



\x\ \x\ 



and 



Hx)-n 00 (^-)\<n 00 (j^-)^-, T>0 n>0; (1.13) 



rr* rf ry 
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In fact and as we shall prove in section 2 this last assumption (jl.l3|) can be largely relaxed if 
for example n — ji m is radial -see (|2.2jh (j2.3J) below. 

We may now state our basic new estimate (already announced in [2B])- Its interest relies of 
course on the bounds stated in Theorem 11.11 

Theorem 1.2 For dimensions d > 2, we assume M.. 9j) . and M.lty and use the notation 

of Theorem M.ll Then the solution to the Helmholtz equation M.l\) satisfies, for R> R and R 
large enough 

12 



M 



[ |V^ 00 (^)| 2 ^^rfx<C Ue+WnlUNnJ^M + M 2 } (1.14) 
J\x\>r Fl Fl L \n i / J 

for some constant C independent of e. 

We would like to point out the sharpness of this inequality. It says that the points where 
I V w n 00 (|^|)| vanishes on the sphere are the concentration directions for the energy \u\ 2 . Indeed, 
we can derive from the Sommerfeld condition below the following proposition. The proof can 
be found at the end of section 3. 

Proposition 1.3 With the assumptions of Theorem \l.J\ below, we have 

liniR^oc— J w^ 2 (| — t)\u{x)\ 2 dx = — 2m j f{x)u{x)dx. (1-15) 

R J\x\<R Fl JR d 



Therefore if u vanishes in a neighborhood of the critical points of we get the bound 

f Hx)\ 2 

I I! dx < oo, 

J\x\>R Fl 

and deduce that Im j Rd f(x) u(x) dx = 0. From the Sommerfeld radiation condition and in 
the constant coefficient case this leads ( [T2j, p. 242) to a restrictive condition on the Fourier 
transform of /, namely /(£) = on the sphere |£| = n 1//2 . It would be interesting to explore 
which should be the conclusion in our setting. A natural hypothesis is to substitute /(£) by 
the generalized Fourier transform defined by S. Agmon, J. Cruz-Sampedro and I. Herbst in 
[2] which depends upon the construction of solutions to the associated eikonal equation. The 
role played by the critical points of was already pointed out by I. Herbst in [TT] . There 
are two related estimates that have been deduced by the same method, the case of stationary 
kinetic equations (the high frequency limit of Helmholtz equations), see [22], and the case of 
two layers, see [1QJ, where the gradient of gives a surface Dirac mass. 

Our next purpose is to give an explicit Sommerfeld radiation condition for the solution ob- 
tained by the limiting absorption principle. Our result complements that of Saito in j25j . 
Therefore we shall assume 

n (x) = X+p(x), A>0, (1.16) 
with p a bounded real function which belongs to C 2 (lR n \{0}) and such that 

\d%p(x)\ < c|a;|- |a| , |a| < 2. (1.17) 
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Theorem 1.4 For dimensions d > 2, assume hl.ify . and $1.16]) . Then for A large 

enough compared to ||p||c 2 ( |x|>i); there exists a unique solution to the Helmholtz equation with 
e = 0, M < oo, and which satisfies for any a > 1, 

[ \Vu(x) - in^uix) ^ | 2 -^r- <C a \ \f{x)\\\ + \x\) a . (1.18) 

jR d \ x \ 1 + \ x \ JR d 

Moreover if there are T > and 5 > such that 

x r 

\n(x) — Wood — r)| < n(x) 1 — tf /or Ixl /arge enough, (1-19) 
\x\ \x\° 

then from il.l(J\) and hi. Iffy we get 

f \Vu{x)-in]i\x)u{x) ^ -^— <C a f |/(x)| 2 (l + Mr. (1.20) 
JR d \ x \ i- ~r \X\ jjgd 



Let us compare the above theorem with previously known results. There is a very extensive 
literature on the limiting absorption principle, see for example 0, [I], [HI, 0> [13, [2Z] and 
references there in. The situation for the Sommerfeld radiation condition is different. When 
n = A + V(x) and V is a short range potential the question was settled by Ikebe and Saito 
in Mochizuku and Uchiyama study in [201 l ar g e range potentials with mild radial oscil- 
lations at infinity like V(x) ~ sin(ln|x|). Hormander in ^H[, chapter XXX, characterizes the 
incoming/outgoing solutions obtained from the limiting absorption principle by some asymp- 
totic behavior, but in his case = A = constant. More general long range potentials were 
considered by Saito in [22] . Although in this latter work perturbations of first order terms 
("magnetic potentials") are also considered let us fix the attention in the conditions for V. 
Saito writes V = p + Q where Q is a short range perturbation, while p satisfies f)1.17jl . Then he 
proves a Sommerfeld radiation condition for A large enough given by Vm ± iv^A(Vy) u, where 

if is an appropriate solution for \x\ > Rq, Rq large enough, of the associated eikonal equation 
i 

|Vvf = l + ^. (1.21) 

Therefore one cannot expect that in general the vector Vy? points at the direction An 
illustrative example is to consider 

p(x) = — —7. 

\x\ 

In this case and for A large enough, see Remark 1.3 in [23], <p{x) — a (^)\ x \ — b(X)xi with 
a(A) = 1/2[(1 + 1/A) 1 / 2 + (1 - 1/A) 1 / 2 ] and 6(A) = 1/2[(1 + 1/A) 1 / 2 - (1 - 1/A) 1 / 2 ]. This 
boundary condition differs from ours in all points except when Vn = (here n = A + p). 
Then the apparent contradiction is clarified thanks to the estimate (|1.14|) which applies for this 
example. 

1 The existence of this solution was established later on by Barles 0]. 
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Notice however that the assumptions for Theorem 11.21 and for Theorem II .41 are different and 
not comparable. In the particular case n = and regular, no smallness assumption is needed 
in Theorem 11.21 because (jl.9|) is trivially fulfilled, while A has to be large to construct the 
solution of the eikonal equation, which is a fundamental step in order to prove Theorem 11.41 
On the other hand Saito's assumption (|1.17j) with n = A +p doesn't imply the existence of the 
limit noo. 

In the proof of Theorem II .41 is crucial the estimate of the tangential part of the gradient given 
in 1)1.10)1 to conclude that for (p given in 1)1.21)1 

/doc 
u{x)\ 2 —rj < +oo. (1.22) 

This is an energy estimate in itself which says that u concentrates along the critical points of 
V T y3. In section 3 we prove that under some conditions these critical points coincide with those 
of Vi-rioo establishing a relation between the energy estimate given in Theorem 11.21 and that 
in (J1.22)) . Notice however that for the proof of Theorem II .21 we do not need the existence of a 
solution to the eikonal equation (jl.21)l which well could not exist. 

The paper is organized as follows. In section 2 we prove Theorem 11.21 Section 3 is devoted 
to the study of the eikonal equation following 4J. In particular we give some properties of the 
corresponding solution. The proof of Theorem 11.41 is given in section 4. In Appendix 1 we 
gather some elementary identities and in Appendix 2 we give the example of the family of wave 
guides mentioned above. 

Acknowledgments. We want to thank T. Hoffmann-Ostenhoff, and G. Barles for enlighten- 
ing conversations. We also thank E. Fouassier for pointing out the shortcomings of a previous 
version. 



2 Proof of Theorem E2 

We are going to prove a more general version. In fact we will consider two different ways of 
measuring |n — rioo|- First recall the assumption (j!.13j) 

\n{x)-n 00 (^ T )\<n{x) T1 , T > 0. (2.1) 
\x\ \x\ 

We can instead assume the alternative conditions 

x r 

\n(x) - n 00 (—)\ < n-—j for Id > -R , T > 0, 5 > 0, and R > 1; (2.2) 

Id \x\° 

there exists (3 < 1, 8 > and T > such that 

(NV T (n-n 00 )-^)_< P\^\ 2 + n(x)j^. (2 ' 3) 
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In particular if n — n M is radial (|2.2|) is sufficient. Also note that from ()1.12|) and either (|2.1j) 
or (|2.2j) we conclude that 

\n\ < C n > — for \x\ large enough. (2.4) 
2 

Let us start recalling the fundamental ingredients of the proof of Theorem 1.1 given in [21]. In 
fact in that paper we discarded the estimate for the tangential derivative because it was of not 
use. However, and as we already said in the introduction, 1)1.11)1 turns out to be fundamental 
in the proof of Theorem 1.2 and Theorem 1.4. 

Consider for R > the special functions <fi given by 

for \x\ < R, 



V* = < \ " for ' i^rV o (2.5) 



(2.6) 



2^r for \x\ < R, 
for | x | > R. 

We will also need the following formulas which hold in the distributional sense 

/' f° r \ x \ — Ri 

i ~ ), > ' ' (^ijlxl 2 - XiXj^j /\x\ 3 for \x\>R, ( "''* ) 



i (d-l)/|x| for |x|>i2. 1 



A simple calculation gives for c£ > 2 



1 vA[2(f>-A^)>-^ vda R . (2.9) 

'5(H) 



4 ,/lR d V r J ~ AR? 



As it is well known A 2 \l/ is positive for d = 2, and has to be treated in a different way. We refer 
to section 5 in [21] for the corresponding modifications. 
Adding the identity (J5.1j) to (|5.3jl with the above choices of ^ and we get 



1 f |Vu| 2 + ! (\Vu\ 2 1 X 



2R.'b(R) J\x\>R^ \ x \ 



d-1 



■ Vm| 2 )/M 



i? 

5(i?) ^ JB(R) 



[R2 , H 2 da R + — I ( //(./■) + ./■ • V//(.r) ) |//| 



+ 1 I A • Vn(x)|M| 2 < 



2 -/|x|>fl l x l 



]- [ f{x)(x-Vu+—±u)+ I f{x)(^--Vu+ d —±u 
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—elm 



(2.10) 



1 f „ - f x ^ -1 

— / X ■ VM U + -r—r- \U U 

■R Jb(R) J\x\>R \ X \ 

We shall estimate separately the various terms in the right-hand side of the above inequality. 
We begin by the four terms containing /. The pairs containing / and Vii is easy to bound by 
the right-hand side of (jl.lUj) because the two norms defined in (jl.4j) and (jl.5j) are one dual of 
the other one (see also the next argument). So we concentrate on the other two terms we have 
for 5 > 

| 7/ | 2 \l/2/ f \ 1/2 

ICQ) Fl 7 V Jc(j) 



< ^ sup 4 / \u\ 2 da R + C 5 N(ff. 

R K JS(R) 



1/2 



(2.11) 



We consider now the terms on n. We have, 



T / (W • Vn)\u\ 2 < V s / raid 



fx • Vn) 



mra 



< ( sup - 



n\u\ 



r> oj +i (ac • Vn) 
> sup 2 J ^ 



< Pi sup — / ;/| // 

B(R) 



x\n r R 



(2.12) 



Finally we have to bound the terms involving e. From the Helmholtz equation one deduces 



\u\ 2 < 



R d 



R c 



\fu\ 



[ |Vd 2 < [ n\u\ 2 + [ \fu\. 

Jn d Jn d Jn d 

With the assumptions ()1.7|) and ()1.8|) we obtain 

/ |Vu| 2 < C( I n 2 \u\ 2 + [ \fu\). 

Using Cauchy-Schwarz inequality and the above inequalities we get 

ef \Vu\ \u\ <Ce 1 f 2 { [ n 2 \u\ 2 + [ \f ul)^ ( [ \f 
jR d VR d M d V VR d ' 

<Ce 1/2 [ \fu\ +C(ssup\n 2 \ [ \f u\ [ W 



1/2 



1/2 



< C(e + sup |n 2 |) 



1/2 



R a 



< c(e + sup i^d 1 / 2 !!!^!!!^-^) 

< SWln^ulW 2 + C s (e + sup \n 2 \) (N(J-)) 2 . (2.13) 

Then plug (|2.11j) - (J2.13)) into (|2.9j) and take the supremum in R. Then we get (ll.lUj) and in 
particular (jl.llj) . 

Let us prove now Theorem 1.2. The proof consists in using the basic equality ()5.3j) with a 
test function that depends on the behavior of n(x) at infinity. We choose for R> Rq such that 
()2.4|) holds and define 

for some non-decreasing smooth function q(r) = for r < 1 and q(r) = r for r > 2. 

With this choice, we will show that the only new information (compared to Theorem II. 1)1 in 
(15. 3|) . is given by the term 

/ Vn(x) ■ W 9 (x) \u(x)\ 2 dx. (2.14) 

JR d 

Then we will take care of the other terms using (jl.K)jh 
As a first step we consider f!2.14|) . We simplify the notation using q = q(^)- Then we get 

f f d \u(x)\ 2 

/ Vn(x) ■ W„(x) \u(x)\ 2 dx = / q I— n^oj)] 2 , cb 
J R d ' J Rd da; 



+ / — n(z) — ^oo(^) |u(x)| 2 <ia;+ / g|x|V r (n(x) -n 00 (— ) ) — n^ta 
J^d or R J R d V |x| / ouo 



l 4^fdx. (2.15) 



The first term on the right-hand side gives the control we look for. As for the second term we 
just have to consider the negative part of J^n(x). Then we proceed as in ()2.12|) to get the lower 
bound 

--^ P ||noo|U« Ill-will/Jo- 

Let us consider first condition (|2.1|l . Then the last term in ()2.15j) is bounded below as follows. 
After integration by parts, it is also given by 

-Tie J (n(x) - n oo(rr)) {plp-oo \u(x)\ 2 + 2 J^ n ooM \x\V T u(x) uj dx 

u(x)\ 2 



> -r 



^oo||c 2 / q n(x) 



\x\ 



dx + 2\\n\\ L o 



q d 

n? \ir~ n oo( u )\ \^ru(x)\ \u\dx 
\xr ouo 



C_ 
R? 



> -— iinooii^nin 1 / 2 ?/!!!! - 



d 2 \u(x)\' 



dx 
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C 2 f \V T u(x)\ 2 
— - n ;« / n-n ax. 



D h--iIL- i | | 

" ■'\x\>R 



Let us assume now ()2.2|) . then the last term of (|2.15j) is bounded below by 

-fif ? |#-nooH| 2 M#^-^|||n^|||| o . 



As a conclusion of this first step we have obtained 

JA) \JL n ( u )\2 \< x> >\ 2 dx 



<Ci/ Vra(x) • W g (x) Kx)| 2 cfe + C 2 M 2 . (2.16) 
The second step is to provide a control on all the remaining terms thanks to the basic equality 

(Q. 

We have 



Vu(x) • D 2 ^ q (x) ■ Wu{x) = — — n^u) |V r w| 2 + -—, n^uj) |— u| 2 

R\x\ R z or 



f Q 

+2 { W\ ~ W ] ne[ ^ u ' VrU] + W VtU ' D " n °° {uj) ' VtM - (2 - 17) 

And because the terms — and q" are supported in the ball {\x\ < R}, we see that from 
(jl.lOj) all the terms in the corresponding integral are bounded by CWriooWc^M/R. 
Next, we consider the term 

A 2 ^ q (x)\u(x)\ 2 = -2Ke / VA^-Vuu 

< 2 / \^A\ Vu \\n 1 l 2 u\. (2.18) 



But we have 



r>l/2 



\VA9 q (x)\ < C Hnoollo-lAf + ^Tial ^ J 172 " ' 



kt: 



3 R\x\ 2 \x\ 2 



Therefore from (Jl.lUj) we get that (|2.18|) is bounded, by 



71, 



di^UcsM/it;. 



n 



The right-hand side terms containing / can be treated as in (j2.11|) and are respectively upper 
bounded by 

CUnoollci^^lllVwIII^, 
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and by 

The last term to be bounded is 



C\\^\\c>N R M)\\\n 1/2 u\\\R - 



elm / Vtf 9 (x) • Vu{x)u{x). (2.19) 
But this is done as in (|2.13|) . The proof of Theorem 11.21 is over. 

3 The eikonal equation 

In order to determine the phase arising in the Sommerfeld radiation condition, we need to solve 
the eikonal equation 

|V<^| 2 = \n(x) = 1 + \p(x), xeR d , A > 0, 

(3.1) 

ip(0) = 0, V(p(x)/\x\ bounded, 

where we assume that A > is large enough and that p 6 C 2 (lR d ) and satisfies 

\D a p(x)\ < C(l + \x\)~ a , for|a|<2. (3.2) 

In order to take into account the linear growth at infinity (which is essential for uniqueness), 
the unique viscosity solution ip to ()3.1|) is better described in terms of the bounded function 
g(x) defined as 

f(x) = \x\g(x). 



This change of unknown yields the following Hamilton- Jacobi equation for g: 
\g\ 2 + 2rgd r g + \x\ 2 \Vg\ 2 = 1 + \p(x), x e R d \{0}, 
9(0) = V^- 

From this H.-J. equation, one can derive properties of the phase (p, which we state now. 



(3.3) 



3.1 Statements of the results 

Let us recall the results of Barles given in [3] . 

Theorem 3.1 (Barles, Under assumption \3. there exists a unique solution to hS.S\l . 
Moreover it satisfies the following estimates, for A large enough and x ^ 0, 

y/l — minp/ A < g(x) < maxp/A, 

\x\\Dg{x)\ < C/X, \x\ 2 \D 2 g(x)\ < C/X, |x| 3 |D 3 ^(a;)| < C/X. (3.4) 
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In the proof of Theorem 11.41 we shall derive the estimate 

\V T (pu\ 2 ^ ^ < +00, (3.5) 

which is similar to our estimate ()1.14|) in Theorem 11.21 Therefore it is a natural question to 
know which is the relation between the two inequalities ()1.14|) and (|3.5jl . In fact we shall see 
that under some extra assumptions on the index n, both inequalities are equivalent. Indeed, 
assume 

\d r p(x)\ < Cr-- 1 " 5 , for some S <E]0, 1[, (3.6) 

with r = 1 2 1 and d r — A • V the radial derivative. Firstly, as a consequence of this assumption, 
n(x) admits a radial limit 

n{r- — r) — > noo(— ) as r — > 00, (3.7) 

and in fact the decay rate in our previous assumption ()1.19j) also follows from ()3.6|) . Then, we 
have 

Theorem 3.2 Under assumptions \3. fy) and \3.b\) . the solution to \3. 3]) satisfies for A large 
enough and 1 ^ the estimate, 

r 1+s \d r g(x)\<C/\, (3.8) 
and, as in \3. 7| ), ^(^jfi) — > ^(jfy) as r —>■ 00, a smooth solution to the equation 

^H^lv^ooMl^nooH, ueS d -\ (3.9) 

Moreover 

\VM = Iv^ooHi + o( r - s ). 



Remark. The solution to the limiting equation ()3.9|) enjoys a specific property. From ()3.4j) we 
get that (#00/ + D^goo) is invertible and after differentiation in ()3.9|) we obtain 

< Ci|V^oo| < (V^riool < CalV^ool- (3.10) 

Therefore the critical points (and thus the extrema) of n M and g^ coincide. 

The derivation of Theorem 13 . 21 can be seen either from the representation through the method 
of characteristics or from the more modern PDE point of view. We give both proofs for the 
sake of completeness. 
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3.2 Proof of Theorem 13.21 based on bicharacteristics 

We decompose the proof in two steps. We first recall the definition of the characteristics and 
how they give a representation formula for the solution to the eikonal equation Ijfl.ljl . then we 
come to the proof of the estimates in Theorem 13.21 

First step. Characteristics. The bicharacteristics are classically defined as the solutions to the 
differential system, parameterized by some q G § d , 



X(t-q) = 


2P(t; q), 


X(t = 


0;g) = 


o, 


P(t;q) = 


±Vn(X(t;g)), 


P(t = 


0;q) = 


itt^vVo), \q\ 2 = l, 


$(t;q) = 


\n(X{t-q)), 


$(t = 


0;q) = 


0. 



(3.11) 



We note for later purpose that 

^[\P\ 2 - \n(X)} = 0, therefore \P(t; q)\ 2 = \n(X{t; q)) Vt G R. 
at A A 

Assume first that we can define a unique diffeomorphism (this involves smallness assumptions) 

x G R d \{0} -> (t,q) G]0,oo[x§ rf_1 , x = X{t;q). (3.12) 

Then, we recall the standard representation formula for the solution (p(x) to (|3.1|) (see [TS] for 
details and complements). 

Lemma 3.3 For allt > and q G S ' -1 such that condition holds, we have 

P(t;q) = V<p(X(t;q)), <p(X(t; q)) = $(t; q), 
for some smooth solution <p(x) to h3. 1)) . 

Proof. Because both the H.-J. solution and the differential system 1)3.11)1 are stable by smooth 
perturbations of n, we can assume that n is constant in a neighborhood of the origin. Then, 
we divide the proof in three steps. 

(i) Then, close to x = the solution to the H.-J. equation is given (p(x) = jj/2\/n(0) \x\ and 
the bicharacteristics are 

X(t; q) = 2t q^n~{ti), P = q^/MO) = V<p(x), $ = 2n(0)t = ^n~{fi)\X\ = <p(x). 

Hence the result holds true. 

(ii) We prove that P(t;q) is a gradient, i.e., is a symmetric matrix. Indeed, from (i), it 
is true for t « and it remains to see that for all times j| Ofjjf) is a symmetric matrix. But we 
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can compute 



A ( 9P) - A 

dt\dx J dt 



ap \ I ax ' 



8(t,q) I \d(t,q) 



dP \ I dX \ 1 , / dP \ d I dX 



d(t,q) J \d(t,q) J 1 \d(t,q) J dt \ d(t,q) _ 

ax w ax \ 1 / ap \ ( ax \ _1 I ax \ i ax 



d(t,q) J \ d(t,q) J \d(t,q) J \ d(t,q) J \d(t,q) J \ d(t,q) 

ap \( ax \ _1 ( dP \( ax N ~' 



\D 2 n{X{t;q)) 

jD 2 n(X(t;q)) - 2. a ^ j j y,; lLin j \ , n! ., n 

iB% ! (.Y(«;,))-2(f)(f) 



because \P\ 2 = \n implies ffjjrj (ff) + ? ' = ix{&\ We nave obtained indeed that 

Jifff) * s a s y mme t r ic matrix and our claim is proved. 

(iii) From step (ii), we can write P{t;q) = Vcp(X(t', g))for some function (p. And since \P\ 2 = 
\n, we obtain that ip solves the H.-J. equation. It remains to identify $ and ip. To do so, we 
write 

j t v{X(t; q)) = V^(X(t; q)) jX{t' q) = 2\P(t; q)\ 2 = jn(X(t; q)). 

This means that ip(X(t;q)) = $(t;g) since the identity holds at t = and the two functions 
have the same derivatives. This completes the proof of Lemma 13.31 

The reciprocal to Lemma f3. 31 is also true. If the H.-J. solution is smooth, then we can derive 
that the bicharacteristic system satisfies the invertibility condition. As before, for simplicity 
we assume that n is constant in a neighborhood of the origin. 

Lemma 3.4 Assume that the solution <p to H.-J. equation is smooth and satisfies x ■ Vp(x) > 
c\x\ and consider the autonomous differential equation 

—X(t;q) = Vcp(X(t;q)), X(0;q) = 0, X(t; q) « n l/2 (0)qt (t small). (3.13) 
dt 

This gives the solution to \3.11\) with the formulas of Lemma \3JA 

Proof. We prove successively that: 

(i) it is the solution to the bicharacteristic system; 

(ii) (t > 0, q E S^ 1 ) i — ► X(t; q) G R d \{0} is one-to-one; 

(iii) it is a diffeomorphism. 

The point (i) follows from the computation in the Lemma f3. 31 and we do not reproduce them 
again. 
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The point (ii) can be proved as follows. 

(a) into : if X(to\ q) = X(t Q + r; </), by uniqueness of the system for backward time X(t; q) = 
X(t + r; q') for all t and thus (considering what happens close to the origin) r = and q = q' . 

(b) onto: take a point y ^ 0, and solve the equation (|3.13j) backward from say t — and call 
Y(t,y) the solution. Because we know that x ■ Vip(x) > c\x\, we have 

dt 2 ~ 1 Wl 

Therefore it reaches K(t; = in finite time (r < 0). Now, the system ()3.13j) is autonomous. 
Therefore we can solve it from t = 0, with q = n(0) -1 ^ 2 lim^o Y(t + s) / s and one has 
V = X(-r;q ). 

For point (iii), we just notice that the above construction, based on the Cauchy-Lipschitz 
theory, also provides C 1 regularity and C 1 regularity of the inverse. 

These two lemmas show that, in the regime of assumption (|3.2j) . with A large enough, the 
solution to the H.-J. equation is given by the bicharacteristics. We can use this fact to conclude 
the proof of Theorem 13.21 



Second step. Estimates in Theorem \3./A We divide the proof of these estimates in two steps. 
First an estimate on the rays is given, then we prove (J3.8|) . Then, the second statement, (|3.9|) . 
is an easy consequence obtained passing to the limit as r oo. 
(i) We have: \X(t; q)\/t = 0(1). Indeed, we can compute 



d \X\ 
~dt~ 

and 



X ■ X = 2X ■ P, 



^-X ■ P = 2\P\ 2 + \x ■ Vn = %n + \x ■ Vn. 
dt AAA 



As a consequence of assumption (j3.fij) . we have 

d 2 \X\ 2 1 



dt 2 4 A 
and thus the result (i) is proved. 



n + o(l) for X large, 



(ii) We have: \x\ 1+s d r g is bounded. Using the calculation of step (i), we have 

q) = I Hn{X(s; q))ds = fifeX ■ P - \X ■ Vn]ds 

= X-P(t-q)-F(t-q) 



with 



F(t; *) = \f X ^ ■ Vn(X( S ; q)) ds, \F(x)\ < C[l + \x\] 



1-6 
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In other words, for x large we have 

and the claim (ii) is proved. This concludes the proof of ()3.8|) and thus the proof of Theorem 
13.21 is complete. 

3.3 Proof of Theorem 13.21 based on H.-J. equations 

The proof is based on the equation for d r g, 

2gd r g + r{d r gf + rgd 2 rr g + r\ Vg\ 2 + r 2 Vg ■ Vd r g = ^d r n. 

Upper bound. We consider the maximum point (if it is not attained, then perturbation methods, 
see HHj, apply) 

max[r 1+5 d r g(x)} = rl +5 d r g(x ), 

and at the point xq, we have 

(1 + 5)^-d r g(x ) + rVd r g(x ) = 0. 

Inserting this in the equation on d r g, we find 

2gd r g + r{d r g) 2 - (1 + 5)gd r g + r\Vg\ 2 - (1 + 5)rVg ■ ^d r g{x Q ) = —d r n. 

\x\ 2A 

Therefore we obtain 

{l-5)[gd r g + T{d r gf\<^d r n : 



and thus we arrive at the upper bound 

max[r 1+S d r g(x)] < ^yy m.ax[r 1+s d r n(x)]/ ming. 



Lower bound. With the same calculation as above, at the point x\ where the minimum is 
attained 

mm[r 1+s d r g(x)} = rl +s d r g(x\), 

we have successively 

(1 - 8)gd r g - 5r{d r g) 2 + r|V^| 2 = ^-d r n, 



2\ 



(1 - 5)gd r g - 5r-^ +28 \r 1+5 d r g) 2 + ^ - g 2 - 2rgd r g] = ^d r n, 

(1 + 8)r 1+s gd r g + 5r-- & {r 1+& d r g) 2 = r 5 [j - g 2 ] - ^d r n, 
which gives the result. 
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4 Proof of Theorem 11.41 and of Proposition 11.31 

We will need the following uniqueness theorem. 

Theorem 4.1 For dimensions d > 2, assume \1. 7[ )-( f7~^) . If u is a solution of Au + nu = 
with u and Vw locally in L 2 , and such that 

liminf / ( | Vw| 2 + \u\ 2 ) do~(x) = 0, as r — > oo, (4.1) 

./[x[=r V / 

t/ien it = 0. 

Proof. Assume -Rj is a sequence going to infinity such that 



linij / (|Vu| 2 + \u\ 2 ) dcr(x) = 0, as j -> oo. (4.2) 



Consider \l/ and as in ()2.5|) and ()2.6|) . Then fix j and use the multiplier V^Vw+(l/2) A^u+<j)u 
in Helmholtz equation 

Au + nu = 0, 

in the region \x\ < Rj, and for R < Rj. Then we repeat the procedure given at the beginning 
of section 2 to prove Theorem 1.1. Note that just the multiplier which involves $ will create 
boundary terms after the integration by parts. Therefore the basic identity ()5.3|) has to be 
modified. In this case we shall obtain 



= / Vu{x)D 2 ^(x)Vu{x)dx 

JB(0,Ri) 



' / \u(x)\ 2 A 2 ^ (x)dx + - / \u(x)\ 2 (Vn ■ V^)(x)dx 



'B(0,Rj) 

'B(0,Rj) * JB(0,Rj) 

+ ~ j d r ^{x)\Vu{x)\ 2 da(x) - K / d r u{x)V^{x) ■ Vu{x)da(x) (43) 
+ ] j d r A^(x)\u(x)\ 2 da(x) - \ I d r (\u(x)\ 2 )A^(x)da(x) 
+ - I d r ^(x)\u(x)\ 2 n(x)da(x). 



2 



Sr- 
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The inner terms are the same as in the proof of Theorem 1.1. Therefore adding those given by 
</> we shall obtain following the same argument given in section 2 



sup — / (| Vu(x)\ 2 dx + n\u(x)\ 2 )dx 

R<R 3 R Jb(Q,R) 

<C f \d r ^(x)\ \Vu(x)\ 2 da(x) + \to I d r u(x)V^(x) ■ Vu(x)da(x) 

R j R 3 



(4.4) 

+ C \d r AV(x)\ \u(x)\ 2 da(x) + C \d r {\u{x)\ 2 )\ \A^(x)\da(x) 

R j R i 



+ C \d r ^(x)\\u(x)\ 2 n(x)da(x). 
Js Rj 

Then taking the limit in j the theorem follows from ([4.2)1 . 

Proof of Theorem 11.41 Let us recall the estimate obtained by Saito in Theorem 1.6 of [25J. 
There is a, with 1 < a < 2 such that 

J \Vu-i\y 2 V V u\ 2 —^^r a <cf \f\ 2 (l + \x\) a dx. (4.5) 

In order to prove this inequality Saito needed the existence of (p, a solution of the eikonal 
equation for \x\ > Rq with Rq large enough, which was established later on by Barles in 0j, 
and that we gave in Theorem 3.1 of section 3. 
From (j4.5|) we get 

J \Vu - a 1/2 V^| 2 ^^ < c J \f\ 2 (l + |x|) a dx. (4.6) 
Hence looking at just the tangential parts of the above inequality and from (jl.ll)) we get 



/ lA^VrH 2 -^"^ / \f\ 2 (l + \x\) a dx. 

J l + \x\ J 



From the eikonal equation we have that 

n-X\dM 2 = |A 1/2 V r¥ f . 
Recall that from the properties of V<£> (|3.4jl we know that d r ip = g(x) + 0(1/ X) > 0. Then 

Also looking at the radial part in ()4.6|) we get 

\d r u - iX l ' 2 d r ipu\ 2 < c j \f\ 2 (l + \x\) a dx. 
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Finally from the above estimates we get 

/ |v u - m 1 / 2 ^u\ 2 = [ (\d r u - in l ' 2 u\ 2 + \V T u\ 2 ) 

y w i + n y 1 + 1 



+ pi 

dx 



< c I (\d r u - i\ 1/2 d rV u\ 2 + \(\ 1/2 d r ip - n 1/2 )u\ 2 + |V r w| 2 ) 



<c I \f\ 2 (l + \x\) a dx. 

Therefore we have proved (|1.18|) . 



Let us assume now that \n — < c(l + \x\) 5 . Then using (jl.lUj) we conclude that 

/ |V« - m^ 2 ^ 2 < c [ \f\ 2 (l + \x\T dx, 

J \x\ 1+ \x\ J 



which is (OHjl . 

Only the uniqueness remains to be proved. From 

\x\ 1 + \x\ 
we get that 

liminf / (\Vu — in 1 ^ 2 —u\ 2 j da(x) = 0, as r —>■ oo, (4.7) 

</|x|=r ^ \ X \ ' 

Multiply in the equation Au + nu = by u and integrate by parts in \x\ < R. Taking the 
imaginary part one gets Im J^ x ^ =r ud r uda(x) = 0. And therefore 

2m u(d r u - in 1/2 u)da(x) = - n 1/2 \u\ 2 da(x). (4.8) 

J\x\=R J\x\=R 

Then from (jl.K)j) . ()4.7j) . ()4.8)1 . and that n is bounded we get that u satisfies ()4.1)1 . Therefore 
uniqueness follows from Theorem 14.11 

Proof of Proposition [T31 The argument is similar to the one to obtain ()4.7j) . We multiply in 

the equation Au + nu = f by u, integrate by parts in the ball of radius Ro and take imaginary 
parts on both sides to obtain 

Im / ufda(x)=Im / u(d r u — in 1 ^ 2 u)da(x) = — / n 1 ^ 2 \u\ 2 da(x). (4.9) 

J\x\=R J\x\=R J\x\=R 

Integrate in Rq in the above expression for < Rq < R, and divide by Rq. Recall that from 
()1.16j) n > c A for some Co > 0. Then the result follows from Sommerfeld radiation condition 
and (fTinjl . 
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5 Appendix 



A.I.- Basic identities. Our proof combines three basic identities that have been used through- 
out this paper and that we state here without proof (see [21] for a proof). For real valued 
functions (p, ip G 5(IR d ), we have 

- f p(x)\Vu(x)\ 2 + - f Ap(x)\u(x)\ 2 + ! ^(x)n(x)\u(x)\ 2 

./rod 2 ./rod ./rod 



TZe / (f(x)u(x)f(x). 



(5.1) 



e I ?p(x)\u(x)\ 2 — Xm / V?p(x) ■ Vu(x)u(x) 



Tm I f(x)u(x)ip(x). 



(5.2) 



Vu(x) ■ D 2 ^(x) ■ Vu(x) - -A 2 ^(x)\u(x)\ 2 + -Vn(x) ■ V^{x)\u{x)\ 2 
-Tie ! f{x)(v^!{x) -Vu{x) + -Am(x)u(x)\ -elm [ W(x) ■ Vu(x)u(x). (5.3) 



A. 2.- Some Examples. In this section we shall give examples of indices n which satisfy 

(x ■ Vn(i))- 



n(x) 



13 < +oo (5.4) 



with (3 as small as wanted and such that Theorem 11.11 does not hold true. In particular the 
estimate for the tangential derivative is false. Also we will exhibit the corresponding Sommerfeld 
radiation condition. 

Condition (|5.4|) is weaker than (|1.9J) and appears naturally in the study of the absence of 
embedded eigenvalues in the continuous spectrum for the Schrodinger operator A + n. Recall 
at this respect the well known example due to Von Neumann and Wigner of a potential which 
satisfies ()5.4j) for (3 large enough and has an embedded eigenvalue, see [21] p. 233. 

We give examples of wave guides which satisfy ()5.4j) but with a scaling which does not leave 
invariant that condition, and therefore there is no possible (3 good for all of them. For these 
examples condition ()1.9|) is not fulfilled either. 

Define Q as the unique positive solution with Q(±oo) = of 



Q" + (q 2 -\)q = 0, ye 
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That is to say Q(y) = seek (y/y/2). Also for A > take Q\(y) = Q(Xy), which solves 

q: + x 2 (qI-1)qx = o. 

Set 6 E C°° a bump function around the origin with 9(x) — if \x\ < 1 and 9{x) = 1 if |x| > 2. 
Then call 

u\{x,y) = Q x {y)9{x)e i ^'\ 
with < e < 1 and Xmy/1 + > 0. Then *u A solves 



G K 2 , 



A< + (n A + iej u\ = f £ (x, y) 
with 

Wa := w A + 1 - -A 2 ; n A := \ 2 Q 2 (\y), 

and 

/ £ (z,y) = hiVTTii sig(x)9'(x)+9"(x))Q x (y)e iVI + r ^. 



Now it is straightforward to check that N(f £ ) < oo. We can pass to the limit in e and we get 
that, setting u x = lim e ^ + u \ — Q\(y)9( x ) e ^ x ^ then 

Au x + n x u x = f(x,y) := (sig(x)9' (x) + f (x^Q x (y)e^. 

Also if < A < 1/2 and = 1 - 4^, we get 

n x -n x 00 =n x = \ 2 Q 2 (\y). 

Therefore 

(&*a) 

r- ^ < 2CA r > 1, 

and can be made as small as wanted. On the other hand straightforward computations prove 
that Hb{r) K| 2 < oo, 

sup J \V(e- iVTTr£lxl u £ x \ 2 < c(l + J {Q'fdy) < c, 

and that there is c > independent of e such that 

/dx 
|V ^ |2 (TT(^TW^- Co|log£| - 

Also notice that ip = iy/1 + ie\x\ is not a solution of |Vy?| 2 = n. 
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